Quasiparticle Structure in Antiferromagnetism around the Vortex and Nuclear 

Magnetic Relaxation Time 
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On the basis of the Bogoliubov-de Gennes theory for the two-dimensional extended Hubbard 
model, the vortex structure in d-wave superconductors is investigated including the contribution of 
the induced incommensurate antiferromagnetism around the vortex core. As the on-site repulsive 
interaction U increases, the spatial structure of charge and spin changes from the antiferromagnetic 
state with checkerboard modulation to that with the stripe modulation. By the effect of the induced 
antiferromagnetic moment, the zero-energy density of states is suppressed, and the vortex core radius 
increases. We also study the effect of the local density of states (LDOS) change on the site-dependent 
nuclear relaxation rate T 1 _1 (r). These results are compared with a variety of experiments performed 
on high T c cuprates. 

PACS numbers: 76.60.Pc, 74.25. Qt, 74.25. Jb, 74.25.Ha 



I. INTRODUCTION 

Much attention has been focused on vortex states in 
various type II superconductors. We are now realizing 
that Fermionic excitations, or low- lying quasi-particles 
induced around a vortex core play a fundamental role 
in determining various physical properties of a supercon- 
ductor and reflect sensitively the pairing symmetryiS*^ 
This is particularly true for <i-wave pairing case realized 
in high T c cuprates. There, the contribution of nodal 
quasi-particles is present in the spectrum of the quasi- 
particles, even before applying a field. 

Recently, several experiments report an "anomalous" 
vortex core in high T c superconductors, which is not ex- 
pected in an ordinary mixed state: (1) Hoffman, et al& 
made an scanning tunneling microscopy (STM) obser- 
vation and found a checkerboard pattern around a vor- 
tex core in Bi22I2, indicating that the low energy LDOS 
is modified with four-site periodicity. (2) Elastic neu- 
tron scattering observes the enhancement of the scatter- 
ing peaks associated with incommensurate antiferromag- 
netism (AFM) modulation in La 2 _ K Sr I Cu0 4 (LSCO)i2^ 
under small fields (~ a few T) compared with H C 2- The 
magnetic moment is found to be increased by a small 
field. (3) Recent /iSR measurement detects a static mo- 
ment in the vortex state in under-doped YBa2Cu30y_,5 
(YBCO)ii£ Kadono, et ali^ also demonstrate that the 
spontaneous moment is induced by a small field around 
a vortex core on LSCO. (4) Finally, site-selective nu- 
clear magnetic resonance (NMR) experiments^ 2 * 1 ^ show 
that nuclear magnetic relaxation time T\ at the core 
site becomes longer as T decreases, indicating the lack 
of LDOS at low T. The above facts strongly sug- 
gest that in high T c superconductors the vortex core 
is quite different from what we expect in "normal" d- 
wave vortex^ namely a picture that under a small ap- 
plied field the spatially modulated AFM moments are 
induced centered around a vortex core and simultane- 
ously the zero-energy peak (ZEP)i4 at a core is sup- 
pressed. This picture is also consistent with the ear- 



lier STM observations which show an "empty core" , 
namely the absence of the ZEP at the core site^ The 
AFM around the vortex 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 was 
also studied theoretically in SO(5) model ) 16 ! 17 QED 3 
models t-J model) 19 i 20 i 21 and Hubbard modeli 23 i 24 i 25 i 26 

Quite recently, two new NMR imaging 
experiments 2 - 7 *^ on the T-dependence of T\ are re- 
ported: Kakuyanagi, et alSL show in Tk^L^CuOg 
that as lowering T, 1/(T\T) at a core site increases 
divergingly toward a temperature Tm below T c and then 
decreases at lower T. This is contrasted with l/(TiT) at 
other sites which exhibit a monotonic decrease below T c . 
This result suggests that below Tm the local magnetism 
appears exclusively at a vortex core site and other sites 
stay in the normally expected d-wave state. Mitrovic, 
et alm& observe a similar divergent behavior in 1/(T\T) 
of YBa2Cu307 at the core site where the crossover 
temperature Tm is quite low, and l/(TiT) at the 
outside sites of the core show a constant at low T. The 
implication of these experiments is twofold: The ZEP 
in d-wave vortex core must be present at Tm < T < T c 
and in some H region, and its contribution enhances 
1/(T\T) in this temperature region. Below Tm, the field 
induced local AFM must exist, which is to remove the 
ZEP, giving rise to a suppression of l/(TiT) at lower 
T. These new experiments further enforce the above 
picture of the AFM-induced vortex core state. 

Here we explore a possibility of the induced incommen- 
surate AFM moment around a vortex in d-wave pairing 
state. It is quite possible because (A) in a d-w&ve vortex 
state the zero-energy quasi-particles produced at a vor- 
tex core are readily available for the magnetism and (B) 
these quasi-particles are under strongly correlated envi- 
ronment in cuprates, where the AFM state appears near 
half filling next to the superconducting phase. Therefore, 
based on these experimental facts we are led to investi- 
gate how two orders; c£-wave superconducting order with 
vortices and the modulated magnetic order can be influ- 
enced in order to compromise mutual competing effects. 

The purposes of this paper are to study detailed prop- 
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erties of the vortex core under the possible magnetic or- 
der, in particular on the quasi-particle structure and to 
calculate the site-selective T\ (r) which turns out to be a 
good probe for examining the quasi-particlcs at a vortex 
core. This study is not only one of our continuous efforts 
on the vortex problem but also associated with our stripe 
problemi22i211 When the stripe structure is formed, the 
doped carriers are accommodated neatly at the stripe re- 
gion, ultimately relating to the mechanism of the high T c 
superconductivity. Therefore, we consider the case when 
the induced AFM around vortex is not in commensurate 
AFM structur e) 25 ! 26 with the ordering vector Q = (tt,tt), 
but in incommensurate spin and charge structure, where 
the spin ordering vector is given by Q = (2tt(^ ~ e),7r) 
or Q = (tt, 27r(i — e)). Following the neutron scatter- 
ing results, 7 ' 8,9 e = | at hole filling ~ |, meaning the 
eight site periodic spin structure. In this case, the charge 
structure characterized by the ordering vector 2Q shows 
four site periodic oscillation, which is consistent with the 
STM observation! 6 . 

The arrangement of this paper is as follows: After in- 
troducing the model Hamiltonian which allows us to de- 
scribe two possible orderings of superconductivity and 
magnetism in an equal footing, namely, the so-called ex- 
tended Hubbard Hamiltonian, we set up the Bogoliubov- 
de Gennes (BdG) equation on a lattice, and a formula for 
the LDOS and Ti(r) in Sec. 2. The basic vortex proper- 
ties, including the amplitude of the induced moment, the 
local charge distribution, and the Fourier component are 
studied in Sec. 3. We discuss the LDOS in Sec. 4, and 
evaluate the T-dependence and site-dependence of Ti(r) 
in Sec. 5. The final section is devoted to conclusion and 
discussion. Some of the results are briefly reported in 
Refs. |23and|H 



II. BOGOLIUBOV-DE GENNES THEORY ON 
EXTENDED HUBBARD MODEL 

We begin with an extended Hubbard model on a two- 
dimensional square lattice, and introduce the mean fields 
n%, a = ( a \ a a i,<?) a * the i-site, where a is a spin index 
and % = (i x ,i y ) and Ag )i)(T = {ai- a a i+& . a ). We as- 
sume a pairing interaction V between nearest-neighbor 
(NN) sites. This type of pairing interaction gives d-wave 
superconductivityi22i22i2i Thus, the mean-field Hamilto- 
nian under H is given by 

i,j.<7 i,<T 
e, 2,(7 

where a\ g (eij )CT ) is a creation (annihilation) operator, and 
i + e represents the NN site (e = ±x, ±y). The transfer 



integral is expressed as 

TT f rj 

t itj = U.j cxp[i— / A(r) • dr], (2) 
00 Jn 

with the vector potential A(r) = ^Hxr in the symmetric 
gauge, and the flux quantum <pQ. The external field is 
introduced as the so-called Peierls phase factor. 

We assume the following hopping integrals iy : For the 
NN pairs tjj = t. For the next-NN pairs situated 

on a diagonal position on the square lattice, » = t'. 
For the third-NN pairs, which are situated along the NN 
bond direction, tij = t". To reproduce the Fermi surface 
topology of cuprates, we set t' = -0.122 and t" = 0.08i. 35 
We consider mainly the pairing interaction V = —2.02. 
The essential results of this paper do not significantly 
depend on the choice of these parameter values. 

In terms of the eigenenergy E a and the wave functions 
U a (r{) and w a ( r i) a t the i-site, the BdG equation is given 
by 

?(°C-^J(":^)^(Sn)- (3) 

where K a ^ t j = — tij + Si j (Urii,— a — /i) with the chemical 
potential /i, Di j — VJ2e A i,j^j,i+e an d a is an index of 
the eigenstatei 32 ! 33 

We study the case of the square vortex lattice where 
the NN vortex is located in the direction of 45° from the 
a-axis, which is known to be a stable vortex configura- 
tion. The unit cell in our calculation is the square area of 
N r 2 sites where two vortices are accommodated. Then, 
the magnetic field is given by H = 2<fio / (aN r ) 2 with the 
lattice constant a. Thus, we denote the field strength by 
N r as HM r - We consider the area of A^^ 2 unit cells. By 
introducing the quasi-momentum of the magnetic Bloch 
state, k = (2tt/ 'aN r N k ){l x ,l v ) : (l x ,l y = 1,2, ••• ,N k ), 
we set u a (r) = u Q (r)e lk ' r , v a (r) = -0 Q (r)e lk r . Then, 
the eigenstate of a is labeled by k and the eigenval- 
ues obtained by solving Eq. within a unit cell. 
The periodic boundary condition is given by the sym- 
metry for the translation R = + lyRy, where 
R° = (aN r ,0) and R° = (0,aN r ) are unit vectors of 
the unit cell for our calculation. Then, the translational 
relation is given by u a (r + R) = u Q (r)e lx ^ r,R ^ 2 and 
v a (r + R) = i a {r)e~ ix{r ^/ 2 with 

X (r,R) = --A(R)-r-2d I (l I -g + -(Hxr ).R 

00 00 

(4) 

in the symmetric gauge. The vortex center is located at 
r + i(3R°+R°). 

The self-consistent conditions for the pair potential 
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Ajj and the number density ni jCT are given by 

Ajj = (ai,;ai,t> = XI u «( r 0<( r j )/(#<*), (5) 

a 

"*,T = <4T a *.t>=Z)M r <)| 3 /(Sa) ) (6) 

a 

™a = <<i«u) = E M r *)l 2 (i - /(£*)) (7) 

a 

with the Fermi distribution function f(E). The charge 
density 7i(rj) = rii,j + rij^, the spin density S z (ri) — 
\{ n i,] ~ n i,l) an d the staggered magnetization M(rj) = 
(— iy x+1, y S z (ri). The <i-wave order parameter at site i is 

A d (r 4 ) = (A Sjl + A_ £:l - Ay,i - A_ flli )/4 (8) 

with A g) j = A i)i+g exp[i^ jfa+ r <+«)/ 2 A( r ) • dr], where 

the singlet pairing component Aj^+g = (oi+g^Oi^) — 
{ai+g^a^i). As the induced order parameter, we can 

R(r,r') 

define the triplet-d-wave order parameter at i-site, given 
by 

A^ lct (r 4 ) = (A* j let + A*_ r g et - A^ lct - A t !^ ct )/4, (9) 

where we use the triplet pairing component A*"^ 1 ?' = 
(ai+i^airf) + {a i+ s^a i: i) instead of A i>i+ g. 

We construct the Green's functions from E a , u Q (r), 
v a (r) defined as 



for the up-spin electron contribution, and 



Ni(E,t) = -Img 22 {r,T,-iu n -^E + i v ) (17) 

7T 

for the down-spin electron contribution. Then, the total 
LDOS is given by 

N(E,r)=N^E,r)+N l (E,r) 

= Y,i\Mr)\ 2 S(E - E a ) + \v a (r)\ 2 S(E + E a )}. (18) 

a 

When we consider the differential tunnel conductance of 
STM experiments, the (5-functions in eq. (|18fl are re- 
placed by the derivative —f'(E) of the Fermi distribu- 
tion function f(E): N(E,r) = - J2 a [K(r)\ 2 f (E a - 
E) + \v a (r)\ 2 f(E a +E)}. _ 

By calculating the spin-spin correlation function 
X+,-(r, r', i£l n ) from the Green's functions in Eqs. i|12|) - 
l|15p. we obtain the nuclear spin relaxation rate^ 

Imx+,- (r, r', iQ n -> Q + i V )/(n/T)\ n ^o 
-EK(r)<(r'K'(r)^(r') 



-v a (v)u* a (r')u a ,(r)v* a ,(r')} 
xirTf'(E a )6(E a ~E a ,). 



(19) 



j(x,x ) 



r\ _ ( 9n(x,x') g 12 {x,x') 



92i(x,x') g 22 (x,x') 



T T [^(x)ip\(x')}) - (T T [ipi(x)ip±(x')]' 



with x = (r, r). After the Fourier transformation of t as 
x') = r^ e --(- T ')j(r, r', Wn ), (11) 

the thermal Green's functions with the Fermionic imag- 
inary frequency u> n = 2irT(n + h) are written as 



(10) 



«a(r)u*(r') 



5n(r,r',cj„) = E 
3i 2 (r,r',w n ) = E 
52i(r,r',w„) = X 

CK 

522 (r,r',c„) = £ 



zw„ — E a 
u a (rK(r') 

v a (r)u* a {r') 

iu> n — E a 

v a (r)v* a (r') 



(12) 
(13) 
(14) 
(15) 



The LDOS is evaluated by using the thermal Green's 
functions as 

Nf(E,r) = —Jmgu(T,T,iu; n -> E + irf) (16) 

7T 



We consider the case r = r' by assuming that the nu- 
clear relaxation occurs at a local site. Then, r-dependent 
relaxation time is given by Ti(r) = l/R(r, r). We use 
S(x) — 7T Im(:r — i?7) 1 to handle the discrete energy 
level of the finite size calculation. We typically use 
i] = O.Olt. In Eq. (|19fl . the first term is proportional 
to iVf(r, E) x Ni(r, E) when r = r'. 

We typically consider the case of a unit cell with 24 x 24 
sites, where two vortices are accommodated. The vortex 
cores are located at (i x ,i y )=(l,l) and (N r /2+l,N r /2+l). 
The spatially averaged hole density is set to = 1 — rTi ~ 
| by tuning the chemical potential fx. By introducing the 
quasimomentum of the magnetic Bloch state, we obtain 
the wave function under the periodic boundary condition 
whose region covers many unit cells. In this paper, tem- 
perature T is scaled by T c which is the superconducting 
critical temperature when {7/4 = 0. 

For the initial state of our calculation, we give a low- 
est Landau level function of the vortex state for A, and 
very weak incommensurate spin modulation for n a with 
the ordering vector Q = (tt, 2tt(| — e)) and e = | appro- 
priate to the hole filling |. By iterating the calculation, 
M{vi) -> when the AFM is absent. When the AFM 
appears, M(r^) grows and makes checkerboard or stripe 
pattern depending on the parameters. 



III. PROPERTIES OF 
ANTIFERROMAGNETISM AROUND THE 
VORTEX CORE 

We solve the BdG equation self-consistently with the 
on-site repulsive interaction [/, using an eight-site peri- 
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odic vertical stripe state as an initial state. We consider 
the case that spin and charge modulations appear around 
the vortex core, which is stabilized by U. Then the self- 
consistent solution of the AFM-induced c£-wave vortex is 
analyzed and compared with the "normal" ci-wave vortex 
for U = O&M 



A. Induced moment around the core 

Figure ^ shows the spatial structure of vortex state, 
i.e., the amplitude of the d-wave oder parameter |Ad(r,-)|, 
the staggered magnetization M(rj), the charge density 
n{vi), the zero-energy LDOS N(E = 0,i"j), within a unit 
cell for various U's. 

First, we consider the induced magnetic moment 
around the vortex core. In d-wave vortex case without 
AFM at U = 0, the superconductivity |A<j(r)| is sup- 
pressed at the vortex core and the zero-energy LDOS has 
peak there (Fig. IHa)). Using these low energy states, 
the staggered moment is induced centered at the vortex 
core for larger U case. It is seen from Fig. [2 that (1) 
At U = no moment appears. (2) The induced mo- 
ment appears centered around the core at U > U cr , and 
the moment grows as U increases. (3) The spatial struc- 
ture of the induced moment changes from a simple AFM 
type to a stripe type when U > 3. It. In the former sim- 
ple AFM state has also checkerboard type modulation. 
This checkerboard modulation is barely seen in Fig. ^c), 
where the "floor" in M(r) is modulated with eight-site 
period. (4) In the stripe case (Fig. ^d)) the amplitude 
maximum of M(r) coincides with the vortex core in this 
example. These tendencies are already seen in our pre- 
vious paper for the vortex state in the stripe stated 

In Fig. Eta), we show the AFM moment at the core 
M v as a function of U for the two cases V/t = — 1 and 
V/t = — 2. The induced moment appears above the criti- 
cal U CT , which depends on V. The critical strength U CI is 
larger for stronger pairing interaction V, as is seen from 
Fig. Et a ); because superconductivity tends to suppress 
the AFM moment. When U <U cr , M v = and the vor- 
tex structure is the same as U = case. We show the In- 
dependence of M v for several U values in Fig. Etb). For 
U/t = 3.4 the moment appears above T c . In this strong 
U case, the spatial structure of the moment is a stripe 
type (see Fig. Etd)), and at T > 0.5T C , the stripe period- 
icity changes from eight site period by the temperature 
effect. For smaller U, M(r) becomes a checkerboard type 
as shown in Fig. [Jc), and M v appears at Tm below T c . 
With decreasing U, Tm decreases. 



B. Charge density around the core 

The spatial variation of the charge density n(r) around 
the vortex core is shown also in Fig. ^ which is also af- 
fected by the induced AFM. While n(r) is suppressed 
at the vortex core at U = 0)2& n(r) increases with ap- 



proaching the core, when M v appears at U > U cx . These 
properties of the charge structure are related to the mod- 
ulation of the LDOS, as discussed later. In the large U 
case giving Tm > T c , the stripe structure becomes emi- 
nent and the charge stripe structure is seen in Fig. ^d). 
When comparing with the spin stripe structure of M(r), 
n(r) has minimum at the node of the M(r) oscillation, 
implying that the excess hole carriers accumulate there. 
The periodicity of the charge modulation is four-site, in 
our doping level ~ 1/8. It is half of the spin structure 
periodicity. 

The {/-dependence of the charge at the core n v is 
shown in Fig. EI c )j where we see that the suppressed 
n v at U < Ucr is increased with the growth of M v when 
U > Ucr- The T-dependence of n v is shown in Fig. Et^)- 
When U = 0, n v is reduced around the vortex core at low 
temperature because of the particle-hole asymmetry^ 
For U > U cr , the decreasing dependence of n v on lower- 
ing temperature changes to the increasing one at T < Tm 
where the moment begins to appear. These charge be- 
havior is similar to the case when commensurate AFM is 
induced around the vortex corei^S& 



C. Fourier component of the spin and charge 
modulation 

In order to see the checkerboard modulation of M(r) 
superimposed in the basic AFM variation, we calculate 
the Fourier transformation S Z q and n q from S z (ri) and 
n(rj), respectively. In addition to the peak of the vortex 
lattice period, S^q has the peaks at q = Qi = (jTt,tt) 
and/or q = Q2 = (w, jtt), where we set the lattice con- 
stant to be unity. These Qi and Q2 are, respectively, 
ordering vectors of the eight-site period spin structure 
along x and y directions. These peaks correspond to the 
neutron scattering peaks observed in LSCOiS*^ As for 
the charge structure, q has the peaks at q = 2Qi and/or 
q = 2Q2. In Fig. Eta), we show the [/-dependence of 
the peak height for S z ,q at q = Qi and q = Q2. When 
U < 3. It, the peak heights at Qi and at Q2 are the 
same, meaning the checkerboard modulation both x and 
y directions. Since the peak height increases, the spin 
modulation becomes eminent with increasing U . When 
U becomes larger than 3. It, one of the peaks for S z ,q 
vanishes, and one of the peaks remains at Qi or Q2, 
which means a one-dimensional stripe, as is seen in Fig. 
dd). 

This checkerboard-stripe transition is also seen in the 
Fourier transformation of the charge structure shown in 
Fig. Et D )- At J7 < 3.1t, n q has the peaks both at 2Qi 
and 2Q2, meaning four site periodic checkerboard charge 
structure. At U > 3. It, n q has the peak either at 2Qi or 
at 2Q2, meaning four site periodic stripe structure. 
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|A d (r)| M(r) n(r) #(r,£~0) 




(b) ?7/t = 2.6 




(d) £//i = 3.4 

FIG. 1: Spatial variation of the vortex state at T = T c /40 for U/t =0.0 (a), 2.6 (b), 3.0 (c) and 3.4 (d). From left to right, we 
show the amplitude of the d-wave oder parameter |Ad(rj)|, the staggered magnetization M(ri), the charge density n(rj), and 
the zero-energy LDOS N(E = 0, rj) within a unit cell. The vortex center is located at the center (i x , i y ) = (N r /2 + 1, Nr/2 + 1) 
and a corner(l,l). 



D. Vortex core radius 

Here, we discuss the vortex core shape and radius from 
the spatial structure of |A<j(r)|. As seen from |A<j(r)| in 
Fig. ma) j when the AFM is absent at U < U CT , we see 
the fourfold symmetric vortex core shape, reflecting the 
fourfold symmetry of the d x 2_ y 2--wSLve pairing potential in 
the momentum spacei After the AFM appears at U > 
U CI , the vortex core becomes circular shape, as shown 
in Fig. dc). In the stripe case at U/t — 3.4, |A<j(r)| 
also shows four-site periodic structure, where |A<j(r)| is 
enhanced at the hole-rich site. In this case, vortex core 
shape is enlarged along the stripe direction. 

To discuss the core radius, in Fig. 0]we show the pro- 
file of the order parameter |Ad(i*j)| and staggered mag- 
netization M(ri) along the ^-direction, which is the next 
NN vortex direction and also the parallel direction to the 
stripe. When the magnetization appears around the core, 
the cJ-wave order parameter is suppressed, therefore the 



radius of the vortex core becomes large. To estimate the 
core radius, we fit the profile of |Ad(rj)| along the next 
NN vortex direction by Am ta.nh A (r / B) , where Am is 
the amplitude in the bulk, and A and B are fitting pa- 
rameters. Our definition of the core radius R is the radius 
where the fitting line is recovered to Am/C = 0.6Am- 
Thus, R is given by R = B In ^{C 1 ^ + 1)/(C 1 / A - 1). 

We show the temperature dependence of the core ra- 
dius in Fig. [S] When U = 0, the core radius R mono- 
tonically shrinks as T lowers, according to the Kramer- 
Pesch (KP) effect^ and saturates at lower T because of 
the quantum limit behavior— even in the d-wave super- 
conductor. For finite U{> U cr ) cases shown in Fig. 
the deviation from the KP linear curve occurs below the 
temperature Tm where M v appears (see Fig. EJb)). That 
is, the shrinkage of R stops at Tm, and the core radius 
increases as T lowers. This is caused by the growth of 
the moment induced around the core which destroys the 
condensation locally, enlarging R. This feature coincides 
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00 



U/t=3A 




TITc 



U/t=3.0 ' U/t=3A 




2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 
f//f 



FIG. 2: The AFM moment M v and the charge density n v 
at the vortex center, (a) [/-dependence of M v at T — T c /40 
for V/t = -1 and -2. (b) T-dependence of M v for U/t = 0, 
2.6, 3.0 and 3.4. (c) [/-dependence of n v at T = T c /40 for 
V/t = -2. (d) T-dependence of n v for U/t = 0, 2.6, 3.0 and 
3.4. 



with the /iSR experiment on L^-xSr^CuC^ by Kadono 
et aZ.fii where as x decreases the R stops at higher T 
and becomes larger. If we interpret the doping depen- 
dence in terms of the [/-change, the smaller x should 
correspond to the larger U. This correspondence is nat- 
ural because as going into underdopings the stripe ten- 
dency is strengthened. Tsuchiura et aZ.j2ij4i assert that 
the core radius becomes smaller as approaching to half- 
filling in their t-J model calculation, it contradicts /xSR 
experimental results and our extended Hubbard model 
calculation. 

The order parameter vanishes at the vortex core 
and recovers to the bulk value continuously. The spatial 
variation of the induces the order parameters with 
other symmetry, such as an extended s-wave or p x ± vp y - 
wave components, around the core^i Here, we only dis- 
cuss the ci-wave triplet pairing component. As increas- 
ing the repulsive interaction U, the magnetization arises 
around the vortex core as mentioned before. In the range 
of U > U c , the checkerboard pattern appears on the floor. 
And in U/t > 3.1, the stripe structure which is modu- 
lated by vortices appears. In these cases, small triplet- 
d-wave component /\*"P lct appears around the core^i^S 
Figure H shows the profile of | A^ riplet (r 4 ) | . By increasing 
U, the maximum point in |/\*"P let | moves to outside the 
core and the amplitudes become large. 



80 



60 



40 



20 



(a) 




; 


c 

- 


•k 


- 




Ac n) 

1 h - 






Stripe 


Checkerboard > 


2.5 3 


3.5 


U/t 





(b) 


I i 


(of) 

- 


- 








(f,0) 


Checkerboard ~>Siripe 

i i i 


2.5 3 


3.5 


U/t 





-20 



FIG. 3: [/-dependence of the Fourier component, (a) Eight- 
site periodic Fourier component of the magnetization, S z ,q for 



Q 



j7r,7r) and Q = (tt, |7r). (b) Four-site periodic Fourier 



component of the charge density, n2Q for 2Q = (^7r, 0) and 
2Q = (0, ivr). T = T c /40. 



IV. LOCAL DENSITY OF STATES 

We study the effect of the induced AFM in the LDOS. 
Figure OJa) shows the LDOS at the vortex core site. 
For U = in the upper panel, we see the ZEP at 
E ~ 0, which is a typical feature of the vortex in d- 
wave superconductors.- The small suppression at E ~ 
may come from the induced other component, such as 
an extended s-wave, around the vortex core4£ However, 
this effect is small. The charge density for each spin is 
given by 

n*,T = / &EN^E,Ti )f(E), 

n u = J dENi(E, Yi)f(E). (20) 

When we see the spectrum inside the superconducting 
gap at U = 0, the weight of N(E, r) for E < is smaller 
than that for E > 0, reflecting the suppression of n(rj) 
around the core shown in Fig. [T^a). 

When the AFM appears at U > U CI , the ZEP is sup- 
pressed around E = Q. With increasing U, the sup- 
pression becomes eminent, and the peaks are shifted to 
larger \E\, as shown in the middle and lower panels in 
Fig. Efa). When M(r<) ^ 0, N^E,^) ^ N^E,^). 
Since n^i < nj.,j in the case S z (ri) < at the site of 
the vortex center, Nf(E,ri) < Ni(E,Ti) for E < and 
Nf(E, Ti) > Ni(E, Ti) for E > 0, following the relation in 
Eq. (J20J. Therefore, split peaks at positive and negative 
energy, respectively, come from the up-spin and down- 
spin contribution. After the AFM appears, the weight of 
the LDOS for E > shifts to that for E < 0, reflecting 
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FIG. 5: Temperature dependence of the core radius R for 
U = 0, 2.6, 2.8 and 3.0. The departure from the monotonic 
decreasing curve of U = 0.0 occurs at the corresponding Tm. 



the enhancement of n(rj) around the core shown in Figs. 
Htc) andHtd). 

In Fig. Etb), we show the LDOS at the site next to 
the vortex center site. There, while the ZEP is smeared 
at U = 0, the suppression of the LDOS around E = 
is seen as in the vortex center site. Since this site has 
up-spin moment, split peaks at positive and negative en- 
ergy, respectively, comes from the down-spin and up-spin 
contribution. The LDOS at farthest site in the midpoint 
between next NN vortices is shown in Fig. EJc), which 
has spectrum as in the uniform <i-wave superconductor 
at zero field. By the induced AFM, the low energy state 
is slightly suppressed. 

The spatial structure of the LDOS at E ~ is also 
displayed in Fig. The ZEP height at the vortex core 
becomes low as the moment M v grows. In Figs. H(c) and 
^d) where the induced moment M v is rather large, the 
peak structure at the core is removed completely and the 
LDOS at the core is small compared with its surround- 
ings, giving rise to a "caldera" type landscape. This is an 



origin of the so-called "empty" core observed by STM, 15 
as discussed later. 



V. NUCLEAR MAGNETIC RELAXATION 
TIME AND THE RESONANCE LINE SHAPE 

As mentioned in the introduction, NMR is a power- 
ful method to simultaneously know the field distribution 
in vortex lattice through the resonance pattern (Redfield 
pattern) and the spatial profile of the zero-energy quasi- 
particles through the relaxation time T\. In this section 
we calculate the site-selective T\ based on the solutions 
obtained in the previous sections and examine the out- 
comes from the "anomalous" ci-wave vortex core. 

In Fig. □ we plot the T-dependence of the 1/{T X T). 
Far from the vortex core (line denoted as "bulk site" in 
the figure), T\ shows the conventional T^ 1 ~ T 3 behavior 
of the line node case as in the bulk at zero field. As 
shown in Fig. Ufa), l/(TiT) shows peak below T c at the 
core site. And, when U/t = 0, 1/(T\T) is increased on 
lowering T due the presence of ZEP. When M v appears 
at T M (<T C ), as shown in Fig. ^b), l/(TiT) behavior 
deviates from the line of the U = case at Tm ■ Below 
Tm, 1/(7jT) is suppressed with lowering T, reflecting 
the decrease of the zero-energy DOS by the splitting of 
ZEP in Fig. Ela). Therefore, we see that, on lowering T, 
1/{T\T) increases at T > Tm and decreases at T < Tm at 
the vortex core site, which is consistent with the observed 
l/(TiT) behavioral At NN-site (Fig. □»), 1/{T X T) is 
also suppressed at T < Tm by losing DOS for E ~ 0. 

Next, we study the internal field dependence of 
1/{T X T). The 1/{T X T) can be distinguished spatially 
by its resonant frequency in NMR experiment. The fre- 
quency depends on the internal magnetic field Hi nt (r), 
which consists of a contribution from the screening cur- 
rent in the vortex lattice and the direct contribution from 
the AFM-moment. At first, we evaluate H? nt (r) from 
the contribution via the screening current distribution 
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FIG. 6: The local density of states N(E,Ti) for U/t = 0.0,2.6,3.0 at T — T c /40. (a) Vortex center site (1,1), (b) Nearest 
neighbor site (1,2) to the vortex center, (c) Bulk site (1,12) in the midpoint between next NN vortices. We also show the up- 
and down-spin contributions, N^(E,Ti) and N±(E,ri). The energy is scaled by t. 
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FIG. 7: Temperature dependence of l/(TiT) for U/t = 0, 
2.6 and 3.0. (a) Vortex center site (1,1), (b) NN-site (1,2) 
next to the center. The l/(TiT) is normalized by the value 
at T c . The line denoted as "bulk site" shows 1/(T\T) at the 
furthest site from the vortex in the midpoint between next 
NN vortices. T M /T C ~0.25 (U/t = 2.6) and 0.75 (U/t = 3.0) 



through the Maxwell equation: 

VxH=— j(r). (21) 
c 

The current is calculated as 

Uri) = 2\e\clm{t i+S:i J2(4+e,* a i,*)} ( 22 ) 

<J 

= 2\e\clm{t i+ e 7 i'^2[u* a (r l+ g)u a (r i )f(E a ) 

a 

+v a (r i+i )v* a (r i )(l-f(E a ))]} (23) 

for the e-direction bond (e = ±x, ±y) at the site r^. 
The screening current jg{vi) circles around each vortex 
core. When M v ^ 0, Jg(rj) exhibits the staggered-like 
fluctuation current^ It does not appear when M v = 
even if U/t > 0. 
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FIG. 8: The distribution of 1/Ti as a function of the internal 
magnetic field H int at T/T c = 0.025, 0.25, 0.5, 0.75 and 0.95 
in the case U/t— 0.0 (a), 2.6 (b), and 3.0 (c). Vertical axis is in 
log-scale, and horizontal axis is in arbitrary unit. The internal 
field distribution function (d) (resonance line shape of NMR) 
at T/T c = 0.025 for 7T lnt (r) (solid lines) and H 3 int (r) (dashes 
lines). The right side (left side) arrow shows the maximum of 

Hint (HU- 



Figure |H1 shows 1/Ti as a function of the internal field 
/Tint for some temperatures. Here, we take account of the 
contribution from the AFM- moment. The internal field 
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is given by 

flint (r) = ^ (/T/ nt (r) + JT ) 2 + m 2 M 2 (r) - 7T 

2 

977 ^ 

" ^« t(r)+ 2i?o M2(r) (24) 

for Hf nt (r) , \mM (r)\ -C i?o- fl) is the external field. 
i/f nt (r) is a contribution of the screening current jg in 
the vortex lattice. The mM(r) with a numerical factor 
to is the direct contribution from AFM moment. Here, 
we assume that the magnetic moments are parallel to the 
CuO-layer as in the AF-state at the half-filling. And we 
choose m to reproduce the observed line shape of NMR. 43 

The internal field distribution function, corresponding 
to the resonance line shape of NMR, is shown for low- 
est temperature case in Figs. OHa)-(c). The horizontal 
width of l/Ti-distribution corresponds to the width of 
the resonance line shape for higher temperature data. 
As increasing temperature, the width of the resonance 
line shape becomes narrow (Fig. Eta)) since the internal 
field narrowly distributes in the vortex lattice. The width 
depends on the radius of the core. As the core radius be- 
comes larger, the maximum point of the screening current 
becomes far from the vortex center, and the peak of the 
internal field at the vortex center becomes broad. Thus, 
the width of the resonance line shape becomes narrow. 

This relation between the width of the resonant line 
shape and the vortex core radius appears also in the U - 
dependence, if we neglect the contribution from M(r). 
We compare the field distribution function of H? nt (r) at 
U/t = 0, 2.6 and 3.0, which are shown by dashed lines in 
Fig. Eld). 111 this figure left side arrows show the maxi- 
mum of H? nt (r). In the case U/t = 2.6 and 3.0 where the 
moment appears around the core at low temperatures, 
the core radius R is larger, and the width of the line 
shape becomes narrower at low temperature. 

The realistic line shape observed in NMR includes the 
contribution of the AFM moment. The solid lines in Fig. 
Efd) show the field distribution function of H lnt (r) given 
by Eq. 1241 In this figure, right side arrows show the max- 
imum of ifi n t(r). In the case U/t = 2.6 and 3.0, by the 
AFM effect, the resonance line shape becomes broader 
and have a long tail toward the higher-TTnt side. After 
the magnetization appears, the ZEP in the LDOS splits 
and 1/Ti becomes small. When we see the 1/Ti-behavior 
as a function of H mt , as in Fig. El 1/Ti shows the same 
behavior as the U = case for T > T M . At T < T M , the 
effective internal field is shifted toward higher- H lnt side, 
and 1/Tj is decreased. These AFM contributions are 
eminent at higher-_f/i nt region, where signals come from 
the spatial region around the vortex core. In Fig. [5Jc) , 
as f/jnt increases, 1/Ti decreases near the vortex core re- 
gion (high H ln t region) at low temperature, reflecting the 
special feature of the spatial structure, namely "caldera" 
structure of N(E = 0,r) shown in Fig. [Jc). This spe- 
cial feature is observed by Kakuyanagi et alM The spa- 



tial structure of this 1/Ti (r)-behavior appears at lowest 
temperature when U is strong (near the stripe state). To 
reproduce the suppression of 1/Ti at the vortex core re- 
gion, the choice of the Fermi surface shape (i.e. t' and t") 
is important. For example, we cannot obtain "caldera" 
type-LDOS structure when t! = 0. 

VI. CONCLUSION AND DISCUSSION 

In this paper, based on the extended Hubbard Hamil- 
tonian, we have studied the vortex structure when the 
incommensurate AFM is induced around the vortex core 
by the competition between AFM and superconductivity 
due to the on-site repulsive interaction U. We calculate 
the spatial distribution of the d-wave superconducting 
order parameter, the magnetization, charge density, the 
local density of states and the nuclear spin relaxation rate 
Ti, and clarify the U- and temperature-dependences. 

As increasing U, the AFM with weak checkerboard 
modulation appears around the vortex core at U > U cr , 
where the AFM- moment appears at Tm(< T c ). At fur- 
ther large U, the checkerboard modulation changes to 
the stripe modulation, where Tm > T c . This transitions 
is clearly seen by analyzing the Fourier component cor- 
responding the ordering vector of the incommensurate 
modulation. 

By the induced AFM around the vortex core, the ZEP 
at the vortex core at U — splits, and the zero-energy 
density of states is suppressed. Reflecting this LDOS 
structure, 1/Ti is suppressed at the vortex core in low 
temperature range. That is, on lowering T, 1/Tj is en- 
hanced at T > Tm and decreased at T < Tm- These 
Ti-behaviors are consistent with the observation^ We 
also discuss the AFM effect on the internal field distri- 
bution, relating to the resonance line shape in NMR. 

The spectrum structure of the LDOS is related to the 
spin and charge structure, such as the suppression (en- 
hancement) of the charge density around the vortex core 
when the AFM is absent (present). The difference of the 
LDOS between up-spin and down-spin is essentially im- 
portant especially to discuss the transport phenomena. 

We also show that the vortex core radius becomes large 
as the AFM appears around the core. This behavior ex- 
plains the doping dependence in recent /iSR experiment. 
In this way, the induced AFM contribution plays an im- 
portant role in the study of the various physical proper- 
ties in the vortex state of high T c -cuprate. 
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